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We present the results of a high resolution UV 2-D spectroscopic survey of 
star forming galaxies observed with HST-STIS. Our main aim was to map the 
Lya profiles to learn about the gas kinematics and its relation with the escape of 
Lya photons and to detect extended Lya emission due to scattering in gaseous 
halos. We have combined our data with previously obtained UV spectroscopy on 
other three star-forming galaxies. We find that the P-Cygni profile is spatially 
extended, smooth and spans several kiloparsecs covering a region much larger 
than the starburst itself. We propose a scenario whereby an expanding super- 
shell is generated by the interaction of the combined stellar winds and supernova 
ejecta from the young starbursts, with an extended low density halo. The variety 
of observed Lya profiles both in our sample and in high redshift starbursts is 
explained as phases in the time evolution of the super-shell expanding into the 
disk and halo of the host galaxy. The observed shapes, widths and velocities are 
in excellent agreement with the super-shell scenario predictions and represent 
a time sequence. We confirm that among the many intrinsic parameters of a 
star forming region that can affect the properties of the observed Lya profiles, 
velocity and density distributions of neutral gas along the line of sight are by far 
the dominant ones, while the amount of dust will determine the intensity of the 
emission line, if any. 

Subject headings: ultraviolet: galaxies — galaxies: starburst — galaxies: halos 
— galaxies: ISM — galaxies: high-redshift 

1. Introduction 

Galaxies with ongoing star formation display characteristic emission lines whose strength 
often dominates the appearance of the optical spectrum (Kcnnicutt, Kobulnicky & Pizagno 
1998). The ionizing radiation from newly formed stars and its interaction with the surround- 
ing gas generate coUisionally excited and recombination lines which become detectable at the 
highest observable redshifts (Melnick, Terlevich & Terlevich 2000). Model spectra of young 
populations predict Lya to be the strongest emission line in the optical/near-infrared (IR) 
spectral range for the simplified assumption of Case B recombination and low metal content 
(see Schaerer (2002) for a very recent set of model calculations). Therefore Partridge & 
Peebles (1967) suggested the Lya line as an important spectral signature in young galaxies 
at high redshift as the expected Lya luminosity could amount to a few percent of the total 
galaxy luminosity. 

Typical Lya fluxes of 10"^^ erg s~^ cm"^ are expected for galaxies at redshifts around 
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3 with star- formation rates of order 10^ yr~^. Such values have been within the reach 
of even relatively modest-sized instruments for several decades. Major observational efforts 
were undertaken to search for Lya emission from such galaxies (Djorgovski & Thompson 
1992). Although quite a few Lya emitters powered by starbursts have been found (e.g., 
Kudritzki et al. (2000), Rhoads et al. (2000), their numbers are generally much lower than 
expected from the observed star-formation rates and Case B recombination conditions. 

The assumption of the Lya intensity as produced by pure recombination in a gaseous 
medium may be too simple. Meier & Terlcvich (1981), Hartmann et al. (1988), Ncufcld 
(1990), and Chariot & Fall (1993) considered the effects of dust on Lya. Lya photons expe- 
rience a large number of resonant scatterings in neutral atomic hydrogen, thereby increasing 
the path length and the hkelihood of dust scattering and absorption. This process can be 
very efficient in removing Lya photons from the line of sight to the observer, leading to much 
lower line strengths in comparison with the idealized Case B. Depending on the aspect angle 
of the galaxy as seen from the observer, this may lead to a decrease of the Lya equivalent 
width. On the other hand, Lya may actually be enhanced due to the presence of many 
supernova remnants which form during the starburst (Shull & Silk 1979). The net result 
is controversial. Bithell (1991) finds supernova remnants to be an important contributor to 
the Lya strength whereas Chariot & Fall (1993) reach the opposite conclusion. 

The theoretical situation is sufficiently complex that observational tests are required. 
The most obvious test arc measurements of Lya in local starburst galaxies whose redshifts 
are sufficiently large to permit observations of their intrinsic Lya outside the geocoronal 
and Galactic interstellar Lya. Observations of local starbursts have indeed been performed 
with the /[/£■ satellite (Meier & Terlevich (1981); Hartmann et al. (1988); Calzetti & Kinney 
(1992); Terlevich et al. (1993); Valls-Gabaud (1993)). Again, the results are controversial. 
For instance, Calzetti & Kinney and Valls-Gabaud find Lya strengths in agreement with 
pure recombination theory whereas Hartmann et al. and Terlevich et al. conclude that 
significant dust trapping of Lya photons must occur. 

The superior spectral and spatial resolution of HST's ultraviolet (UV) spectrographs has 
allowed new insight into the formation process of Lya. HST-GHRS spectroscopy of eight 
gas-rich irregular galaxies by Kunth et al. (1998) indicates yet another, and most likely the 
dominant parameter governing Lya emission: neutral gas kinematics. Kunth et al. found 
Lya emission with blueshifted absorption in four of the galaxies. In these objects the 01 
and Sill absorption lines are also blueshifted, suggesting an outflow of the neutral gas with 
velocities around 200 km s~^. The other four galaxies show broad damped Lya absorption 
profiles centered on the wavelength of the ionized gas with no detection of Lya emission. 
The observed galaxies span a metallicity range of more than a factor of 10 and display no 
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correlation between metal abundance and Lya emission strength, a correlation that had 
been postulated on the basis of the lUE spectra and on theoretical grounds due to the 
appearance of dust. The velocity structure of the neutral gas in these galaxies is the driving 
factor that determines the dctectability of Lya in emission. When most of the neutral gas 
is velocity-shifted relative to the ionized regions, the Lya photons can escape, a suggestion 
supported by recent models of Tenorio-Tagle et al. (1999). Nevertheless, while the velocity 
structure of the neutral gas is the driver for the detectability of the Lya emission, we want 
to stress that the amount of dust will still be responsible for the intensity of the line, when 
observed. The Lya photons not affected by resonant scattering will still be strongly affected 
by dust extinction, which has a maximum in the UV range, around the Lya wavelength. 
The properties of the neutral gas (density, kinematics, covering factor) determine thus the 
shape of the profile and the equivalent width of the line, while the amount of dust drives 
only its intensity. The implication is that feedback from the massive stars via ionization 
and the creation of supcrbubblcs and galactic-scale outflows lead to the large variety of Lya 
profiles. The escape of Lya photons depends critically on the column density of the neutral 
gas and dust, the morphology of the supershells, and the kinematics of the medium. Since 
these effects can be highly stochastic, theoretical predictions for the Lya strength are quite 
uncertain, and empirical guidelines are called for. 

A similar pattern seems to apply to high redshift galaxies. At 2.5 < z < 5.2, about half 
of the galaxies found show Lya in emission (Steidel et al. (1999); Rhoads et al. (2000); Prye, 
Broadhurst & Benitez (2002)). At the same time, the Lya line is asymmetric, displaying 
an extension towards larger wavelengths. The blue "edge" of the fine could be described as 
showing a P-Cygni profile, and the centroid of the line is redshifted by some hundred km s~^ 
with respect to the metallic absorption lines. This is consistent with gas outfiows, breakout 
of the gas bubble produced by the star forming region, and results in the escape of Lya 
emission. 

In order to analyze the spatial (kinematical) structure of the Lya emission in nearby 
star forming galaxies, we used the Space Telescope Imaging Spectrograph (STIS) on board 
HST to reobserve three galaxies of the Kunth et al. (1998) sample, two of which clearly 
showed with GHRS a P-Cygni Lya profile (Haro 2 and IRAS 0833-1-6517) and a third one 
that only presented a broad Lya absorption (IZwlS). We aimed also to analyze whether Lya 
could be leaking in external regions, while being completely absorbed around the starburst 
itself. 

The data and data analysis technique are presented in Section 2; the expected structure 
of an HII region and subsequent effects on the resulting Lya profiles follow in Section 3; the 
superbubble generated by the HII region and their time evolution are analyzed in Section 4, 
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together with the expected effect on the observed Lya hne profile. The general discussion is 
presented in Section 5, and the particular implications for high redshift galaxies, in Section 
6. Finally, Section 7 of the paper presents the summary and conclusions of our work. 

2. HST-STIS observations 

We list in Table 1 the basic data for the galaxies observed with HST-STIS. The journal 
of observations is summarized in Table 2. We will complement our discussion with some data 

presented by Thuan & Izotov (1997). These autors analyzed a sample of 3 compact galaxies 
experiencing strong starbursts. They found Lya emission in only one case, Tol 1214-277. 
The line was very strong and rather symmetric with no evidence of blueshifted absorption. 
They concluded that in this case the line was visible because the covering factor of the 
neutral cloud was small, leaving Lya photons able to escape through paths relatively free of 
HI gas. 

Previous GHRS spectra were taken through the Large Science Aperture (LSA), com- 
prising roughly 1'.'7 x 1'.'7 (2"x2"for the pre-Costar IZwl8 data), centered on the maximum 
UV continuum. STIS data have been obtained through the 52"x0'.'5 slit on the FUV MAMA 
and optical CCD detectors. 

In the FUV range the G140M grating was used. The data have been extracted with 
standard STSDAS procedures. The different integrations obtained were combined and av- 
eraged with the MSCOMBINE task to improve the signal to noise ratio. The detector and 
geocoronal background was removed averaging 100 rows on areas below and above the stellar 
continuum. While the detector background was neghgible for Haro 2 and IRAS 0833-1-6517, 
it affected significantly the bluest continuum of IZwlS. 

Low resolution optical spectra were taken with the STIS CCD using the G430L grating, 
in order to compare locally the extension and intensity of Lya and optical emission lines, 
namely [OH], [OIII] and Hp. The optical spectrum of IZwlS was underexposed, but we 
could clearly detect the extension of the optical lines for Haro 2 and IRAS 0833-1-6517. 

The UV spectral images have a spatial scale of 0'.'029 and a spectral dispersion of 0.053 
A/pix, giving a spectral resolution around 0.15 A, translating to ~ 37 km s~^ at the Lya 
wavelength. The optical images have correspondingly scales of 0'/05 and spectral dispersion of 
2.7 A/pix (spectral resolution ~ 6.7 A). We have multiplied the resulting fluxes by the STIS 
parameter diff2pt (diffuse to point source conversion factor for absolute photometry), as 
required for non-extended sources - although resolved, the objects are indeed rather compact. 
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In the next section we describe in detail the observations object by object. 

2.1. Haro 2 

We show in Fig. 1 the STIS spectral image of Haro 2 with different scales. The slit is 
oriented along the minor axis of the galaxy (during a first observation with the slit along the 
major axis the aperture was unfortunately misplaced due to an observer error). It can be 
seen that the UV continuum is very compact, extending over only around ff!9 (corresponding 
to ~ 85 pc). There is a very strong, quite extended and spatially asymmetric Lya emission 
line. Lya emission is detected over 7", 6" of which are located to the upper (SW) part of 
the slit. There are hints that the nebulosity is even more extended, at a fainter level, to this 
direction. The peak of the Lya emission is slightly offset (by 0'.'2 or 19 pc) from the peak of 
the UV continuum. A spatial profile of the Lya emission and the UV continuum is shown 
in Fig. 2. In Fig. 3 we display the extracted spectrum of the central region of this galaxy. 

The most striking result from these data is that the conspicuous and completely black 
absorption edge affects the blue wing of the Lya profile over the whole region where Lya 
is detected, at essentially the same velocity. This implies that the neutral gas that causes 
the resonant scattering is approaching us at around 200 km s~^ and extends over at least 
~8"(~750 pc) with no hint of any velocity structure. We show in Fig 4 the Lya profile 
at different positions over the slit. The zero in the velocity scale corresponds to the Lya 
centered at the redshift of the central HII region determined from the optical emission lines. 
These plots confirm that the blue edge of the Lya profile always appears at basically the 
same velocity, with variations smaller than ± 30 km s^^. Note as well that the red wing of 
the Lya emission profile extends essentially to the same velocity (~ +500 km s~^) all along 
the slit, independently of the peak intensity of the line. As we will discuss later the emission 
from this wing originates from a receding shell far from the central HII region. 

The optical spectral image is shown in Fig. 5. The optical continuum, from 3000 to 5000 
Ais very blue and appears to be dominated by very young stars, as illustrated in Fig. 6. The 
already known Wolf-Rayet feature at around 4686 A is clearly detected in the spectrum. 
Emission lines are detected over ~l'/5 extent, which corresponds to the core of the Lya 
line. The optical lines show some asymmetry at the upper part of the slit. The spectrum 
has enough signal-to-noise so as to measure the total flux of HP and the OII+OIII lines; 
these are listed in Table 3, together with the Lya fiux. It is important to note that these 
fiuxes are integrated over the same area (l'.'5x0'.'5) hence their ratios refiect the intrinsic 
local ratios in the gas. The observed Lya/ Hj3 ratio is only ~2.1, much below the theoretical 
value Lya/Hl3r^33. Mas-Hesse & Kunth (1999) measured an extinction of E{B — V)=0.22, 
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as derived from the Balmer lines observed through a rather large aperture containing the 
whole ionized region. With this extinction, the expected Lya/H(3 ratio would be around 
3.4 (assuming a Large Magellanic Cloud (LMC) extinction law in the UV - see Mas-Hesse 
& Kunth (1999)). This implies that the effect of resonant scattering has lowered the Lya 
intensity by at least 40%. Nevertheless, we stress that this extinction value represents an 
average over a large region. Since the Balmer line ratio is not measured with the same 
spectral resolution as the Lya line, it is not possible to disentangle effects due to dust 
extinction with those of resonant scattering. 

This galaxy exemplifies the paradigm defining the dctcctability of the Lya emission 
line: while most of the line emission has been destroyed in this case by dust absorption, the 
line is still detectable, and indeed quite prominently, due to the kinematical configuration 
of the neutral gas surrounding the HII region. If there were no dust at all, the line would 
remain undetectable unless the neutral gas were expanding. On the other hand, if there 
were not neutral gas at all, both the Lya emission line and the surrounding continuum 
would experience roughly the same absorption, so that the equivalent width of the line would 
remain unaffected, independently of the amount of dust. Only the absolute line intensity 
would be strongly dependent on the presence of dust. 

The evolutionary synthesis modelling by Mas-Hesse & Kunth (1999) yielded an age 
for the starburst episode around 4.8 Myr and more than 6.6 x 10^ M0 of gas having been 
transformed into stars. As cited by these authors, there is evidence in Haro 2 of previous 
recent starburst episodes, which could have taken place in the last 50 Myr. 

2.2. IRAS 0833-1-6517 

The STIS spectral images of IRAS 0833-1-6517 reveal a much more complex structure 
than Haro 2 (see Fig. 7). The UV continuum shows a patchy distribution, and is extended 
over ~3" (1.2 kpc) - its limits arc diffuse: the UV continuum seems to be even more extended 
at a lower brightness level. There is also a strong, quite extended and spatially asymmetric 
Lya emission line. Lya emission is detected over at least ~10" (4.1 kpc). Lya emission 
shows 2 main peaks: the brightest is shghtly offset with respect to the center of the UV 
emission, ~0'/3 (124 pc) to the upper (S) part of the slit. The second peak is also offset, 
~1'.'2 (496 pc) to the lower (N) part of the UV continuum. Some weaker peaks can also be 
identified on the image. This is shown in Fig. 8, where we plot the spatial profile of the 
Lya emission, on top of the continuum profile. The extracted spectrum in the Lya. region 
is given in Fig. 9, showing the asymmetric Lya emission line with a broad absorption to the 
blue. 
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We show in Fig. 10 the Lya emission profile across different regions. As already found 
in Haro 2, the blue absorption edge has no any significant velocity structure over at least 
the 10" (4.1 kpc) along the slit where emission is detected (the velocity at which the profiles 
go to zero in the different regions is within the range —50 — 1-50 km s^^. Here again, the 
neutral gas in front of the Hll region, within at leat 4.1 kpc, is approaching us at basically 
the same velocity, ~300 km s~^. Moreover, the red wing of the emission profile extends to 
roughly the same velocity, around 700 km s~^, as also found on Haro 2. 

We confirm the detection made on the GHRS spectrum of a secondary Lya emission 
line, blucshiftcd by around —300 km s~^ with respect to the Hll region velocity determined 
from the optical emission lines. This secondary emission shows some interesting features: 

• It is clearly offset from the continuum, peaking at ~0'.'9 (370 pc) from the UV con- 
tinuum, to the upper (S) part of the slit. The emission, at a weaker level, extends 
clearly to the lower (N) part of the sht, at least over the region where there is some 
UV continuum. 

• Although the signal to noise is low, this secondary emission peak shows some velocity 
structure, with components rcdshifted by ~100 km s^^ from its average centroid. In 
general, the profile is not a gaussian one and seems to be the convolution of various 
emissions at different velocities. 

We show in Fig. 5 the spectral image of IRAS 0833+6517 taken in the optical range. 
The continuum has a narrow peak about ~0^'3 wide (124 pc), on top of a weaker extension 
over ~2"(826 pc). There is a region to the upper part of the sht which is clearly more 
conspicuous in the UV than in the optical. It is here where Lya and the forbidden oxygen 
lines peak. This region is offset by ~1'.'2 (496 pc) from the optical continuum maximum. We 
show in Fig. 11 the optical spectra of the bright nucleus and of this region. 

We believe that the maximum of the optical continuum is due to relatively older stars, 
with a weak contribution to the UV and to the ionization. Indeed, coincident with the optical 
continuum maximum there is a local minimum in all emission lines strengths. The optical 
continuum shows a very prominent Balmer decrement, as well as Balmer absorption lines 
of stellar origin, but only in the central region (see Fig. 11). The optical continuum in the 
upper part, where the emission lines peak, is much bluer and with weaker Balmer decrement 
and Balmer stellar absorptions, as shown on the figure. Gonzalez Delgado et al. (1998) 
presented a detailed analysis of the massive young stellar population in this galaxy based 
on evolutionary synthesis models. They derived a relatively old age for the burst between 
around 6 and 7 Myr. These authors concluded that there was a significant dilution by an 
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underlying older stellar population, which we identify indeed with the population producing 
the maximum of the optical continuum. 

Finally, [Oil], [OIII] and if/3 are detected over ~3'/5 (1.4 kpc) and coincide again with 
the core of the Lya emission. The measured Lya/H(5 ratio is ~2.1, indeed similar to the 
value measured in Haro 2 (see Table 3). Gonzalez Delgado et al. (1998) quote a large range 
of extinction values for this galaxy, with E[B — V) between 0.17 and 0.52. According to these 
authors, the value derived from the Balmer lines ratio is E{B — V)= 0.52. Assuming this 
extinction and an LMC law, the expected Lya/Hj3 ratio would be only around 0.28, much 
smaller than the value measured by us. On the other hand, if wc consider the extinction 
value derived from the continuum fits, with E[B — V) = 0.17, and the same LMC law, the 
expected ratio would be much higher, Lya /HP ~ 7, indeed significantly higher than the 
measured value. We conclude again that although the effects of extinction and scattering 
can not be disentangled, our results show clearly that both effects must significantly affect 
the observed properties of the Lya emission line. 



2.3. IZwl8 

We show in Fig. 12 the STIS UV spectral image of IZwl8. The UV continuum is quite 
extended (over ~2", corresponding to ~97 pc). The total flux within the STIS aperture is 
around a factor 4 lower than within the GHRS aperture (2"x2", pre-Costar data), which is 
consistent considering the extension of the source. 

A very broad and damped Lya absorption is detected, with complete blackening at the 
center of the line, as shown in Fig 13. The longer wavelength range of STIS allows for a 
more complete coverage of the absorption red wing, as compared to the GHRS data from 
Kunth et al. (1998), confirming the properties of the damped absorption profile reported 
there. No emission is seen at any position along the slit and regardless the intensity of the 
continuum. This is clearly understood if the neutral gas covers the whole ionized region 
with a large column density. This neutral gas must be static with respect to the central 
HII region all along the region covered by the slit, so that resonant scattering affects all 
Lya photon emitted by the central HII region. The complete blackening of the absorption 
profile indicates that essentially all photons have been finally destroyed, most likely by dust 
absorption. 



-10- 



3. The structure of an HII region and effects on the resulting Lya profiles 

An HII region surrounding a cluster of massive newly formed stars should produce 
an intrinsic Lya emission line with a Lya/ Ha ratio around 11 and an intrinsic velocity 
dispersion (a) similar to the optical Balmer lines one. Nevertheless, there are several factors 
affecting significantly both the intensity and the profile of the Lya line as observed from 
outside the HII region. We will discuss the different effects in this section aiming to explain 
the observed Lya profiles. 

The main effects contributing to modify the Lya emission line are the following: 

• Absorption by dust. The known extinction curves peak in the far UV range (around 
1000 A), so that the interstellar extinction will be maximum around Lya. As an 
example, an E{B — V) of 0.2 will yield the absorption of 38% Ha luminosity, but of 
96% of the Lya emission - according to the Small Magellanic Cloud law, generally 
valid in these environments (Mas- Hesse & Kunth 1999). 

As a first approximation, the Lya line and the surrounding continuum should experi- 
ence the same extinction by dust (assuming there is no neutral gas along the line of 
sight), so that the equivalent width of the line would not be affected. Therefore, in 
regions where the continuum is well detected, dust absorption alone could not explain 
the weakness or even absence of Lya in emission. 

Nevertheless, we want to remind that the spatial distribution of stellar continuum 
sources, ionized gas and dust clouds might be very different, as found in nearby star- 
burst galaxies (see for example the analysis of NGC 4214 by Mafz et al. (1998)). This 
spatial decoupling might lead to very different extinctions affecting the UV continuum 
and the emission lines, as discussed by Mas- Hesse & Kunth (1999). 

• Scattering by neutral hydrogen. Lya photons travelling through neutral hydrogen 
will suffer resonant scattering which will re-distribute them within the cloud. This 
effect dramatically increases the sensitivity of Lya photons to dust absorption. In 
the presence of neutral hydrogen we therefore expect the complete destruction of Lya 
photons by dust absorption, even in environments with relatively low dust abundance. 

• Presence of expanding shells. The presence of an expanding shell around the HII 
region can dramatically affect the shape and intensity of the Lya emission line. We 
have identified 4 major elements: 

— If the expanding shell is sweeping the neutral gas, the resonant scattering will 
affect photons with energy slightly higher than those of the central Lya emission 
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line. As a result, only part of the Lya emission will be able to go through the 
neutral medium and become visible. P-Cygni like profiles will be expected in 
these cases, as we will show later. 

— If a fraction of ionizing photons escapes the central HII region, the ionization 
front could reach the internal layers of the shell. Wc would expect to detect 
then a redshifted Lya emission component originated at the inner surface of the 
receding shell. 

— In addition, part of the Lya photons produced in the central HII region could be 
backscattcred by the neutral layers of the receding shell to the obsever's line of 
sight, appearing also redshifted (with respect to the HII region systemic velocity) 
by the shell expansion velocity. 

— Finally, a recombining region, behind the leading shock, could cause a secondary 
Lya emission component, which would be observed blueshifted by the shell ex- 
pansion velocity. Since this component originates outside the neutral shell, it 
wouldn't be affected by resonant scattering by the neutral layers. 

We will analyze in the next sections in more details effects related to resonant scattering 
and the presence of an expanding shell. 

3.1. Effect caused by resonant scattering 

Neutral hydrogen will scatter resonantly the Lya photons inside the cloud. In Fig. 14 we 
show the expected absorption profiles produced by a layer of neutral hydrogen with different 
column densities. Profiles have been computed using the XVoigt code (Mar & Bailey 1995), 
assuming a thermal broadening of the neutral cloud with 6 = 20 km s~^. The Figure 
shows how huge are the effects of resonant scattering by neutral hydrogen: for relatively 
small column densities (A^ lO^'^ cm'^) the center of the line reaches total absorption. 
For column densities higher than N ~ 10^^ cm~^, the line becomes damped and the full 
saturation (hence no transmission) spreads rapidly towards both sides. Remark from the 
Figure that for A" ~ 10^^ cm~^ all photons emitted by particles within a range in velocity 
of ± 1000 km will be completely absorbed. This not only affects the Lya emission line 
photons, but also those photons from the surrounding continuum emitted by the central 
cluster of massive stars. Moreover, the wings of the absorption would affect photons at 
energies corresponding to ± 4000 km s~^ from the line center. 

In principle resonant scattering by neutral hydrogen does not destroy the incident Lya 
photons. In a completely dust-free cloud, photons would be internally scattered until they 
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reach the external surface of the cloud from where they are able to escape. We would expect 
Lya photons to leak from the overall surface of the neutral cloud surrounding the HII region, 
producing a very low surface brightness. Ahn, Lee, & Lee (2001, 2002) have discussed the 
effects of Lya photons scattering by neutral Hydrogen in a dustless and static medium. 
They assume that no photon become destroyed (in absence of dust) and then analyze the 
redistribution in frequencies of the escaping radiation. They predict the formation of a 
rather narrow absorption trough at the line center, with extended red and blue wings but 
their profiles differ significantly from the observed ones basically because the effect of dust 
was not considered in these works. Indeed Kunth et al. (1998) showed that when a static 
cloud of neutral gas is surrounding the HII region the resulting profile corresponds to a typical 
Voigt-like absorption line. We stress again that even very small amounts of dust completely 
destroy all scattered photons via multiple scattering events. Our simulations consider that 
all photons affected by scattering, according to the corresponding Voigt function, will be 
finally absorbed by dust and will not be re-emitted. 



3.2. Effects caused by an expanding shell 

According to the profiles shown in Fig. 14 it would become nearly impossible to detect 
Lya in emission from starburst galaxies, since they are usually immersed in rather dense 
neutral clouds. However, we argue that the presence of an expanding shell, with properties 
evolving concurrently with the central cluster of massive stars, allows for an explanation of 
the various profiles actually detected in most objects, as discussed in Tenorio-Tagle et al. 
(1999). In the following discussion the expansion velocity of the shell, v^xp-, will refer to the 
velocity of the shell with respect to the central static HII region. 

We first show in Fig. 15 the expected Lya profiles when the neutral gas surrounding 
the HII region is expanding at a certain velocity {vexp)-, for different HI column densities 
along the line of sight. As explained above, resonant scattering in the expanding neutral 
gas affects photons with higher energies than those emitted in the Lya line by the central 
HII region. As a result, the absorption of Lya photons mostly affects the blue wing of the 
emission profile producing a classical P-Cygni shape. The damped part of the absorption 
profile remains completely black and the emission line becomes strongly asymmetric. If the 
HI column density is high enough (in our example, above around 10^^ cm^^), the absorption 
is total and only a broad damped profile would be detected. The resulting spectrum will 
nevertheless depend on the convolution of different parameters: expanding velocity of the 
shell, HI column density along the line of sight, intrinsic intensity and width of the central 
Lya emission fine, etc... so that we expect to detect the different cases presented in Fig. 15 
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almost randomly. 

Note that P-Cygni Lya emission profiles show a peak which appears redshifted with 
respect to the intrinsic line by up to several hundred km s~^. This is shown as an example 
in Fig. 16. The broader the intrinsic Lya emission line, the redder the peak can appear after 
partial absorption by an expanding shell. Clearly, the effect of resonant scattering within 
a dense, expanding neutral shell leading to a redshifted Lya emission peak, should not be 
taken as evidence for the presence of outflowing ionized gas. 

In the presence of a neutral expanding shell we expect to detect interstellar absorption 
lines associated to the neutral or weakly ionized gas. These will appear blueshifted with 
respect to the systemic velocity by —Vexp, as indeed found by Kunth et al. (1998). 

If the expanding shell remains symmetric, the internal surface of the receding section 
of the shell will also affect the observed Lya profile: a fraction of the ionizing photons can 
ionize this internal surface of the shell. The photons produced by the approaching section 
of the shell would be immediately scattered by the neutral gas. Since both the ionized and 
neutral layers are expected to expand at the same velocity, the absorption would be very 
efficient and no photon can escape. On the other hand photons originated at the receding 
section of the shell would be emitted with a redshift corresponding to +f eaj^and will be able 
to travel freely through the approaching neutral layers of the expanding shell. 

Moreover, a fraction of the photons scattered by the neutral layers on the receding shell 
will end being rc-cmittcd on the direction of the line of sight, also with a redshift of +Vexp, 
and able to go through the neutral layers of the approaching shell. As a result, we expect 
to detect an additional Lya emission component peaking at +Vexp- This component should 
be broader than the intrinsic Lya line, since the resonant scattering affects photons over a 
rather wide range of energies. Nevertheless, we expect this component to be rather weak 
(at a level of few percent of the intrinsic Lya line) producing just a broadening at the red 
wing of the resulting Lya emission profile since most of the photons affected by scattering 
would be destroyed by dust. This effect is shown in Fig. 17, where we have plotted the 
expected resulting profile assuming a component originated at the receding shell (both by 
photoionization and backscattering) amounting to 10% of the intrinsic Lya emission line. 

Ahn, Lee, & Lee (2003) have extended their models Ahn, Lee, & Lee (2001, 2002) to 
a spherical, expanding supershell. As discussed in Sect. 3.1, these models do not consider 
the effects of dust absorption, and only allow for multiple backscattering events on the 
internal layers of the expanding neutral shell, both on the approaching and on the receding 
sides. After each backscattering photons get an additional redshift of +Vexp but since no 
photon destruction is considered several redshifted emission peaks are produced by multiple 
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backscattering. However we argue that when absorption by dust is taken into account, 
multiple peaks anticipated by these authors would disappear. Moreover, the intensity of 
the single rcdshiftcd peak would be much weaker than their prediction since most photons 
affected by scattering are destroyed before being back-emitted. We have not performed a line 
transfer model to reproduce this effect. We just qualitatively show that both backscattered 
photons and the Lya ones originating at the internal, ionized layers of the expanding shell, 
will produce an additional weak and broad emission component centered at around +Vexp- 

Finally, an ionized region can develop at the external shock front of the shell. The 
recombining medium immediately behind the shock front would produce an additional Lya. 
emission line unaffected by resonant scattering. This emission line observed at a blucsliift 
corresponding to the expansion velocity of the shell, —v^xpi if present, would fill partially the 
damped absorption profile. Its redshifted counterpart, on the other hand, will be scattered 
and fully absorbed within the large column density of neutral matter in the receding part of 
the shell. 

The above discussion aims to confront the reader with the view that there are many 
parameters related to the structure of the HIT region that affect the resulting Lya profile. In 
particular, the velocity and density structure of the neutral gas along the line of sight seems 
to be a dominant factor for the visibility of the line. 



4. Lya emission in starburst galaxies: an evolutionary view 

In realistic cases, bubbles and superbubbles generated around an HII region will evolve 
with time, in parallel with the evolution of the massive stars themselves. As a result, the 
physical conditions that the Lya photons experience when escaping from the HII region will 
differ drastically. Therefore the expected Lya. profiles will be strongly dependent on the 
geometry and evolutionary state of the starburst. 

We show in Fig. 18 the predicted evolution of the expanding shells generated by an HII 
region as they interact with the disk and the halo of the host galaxy. The figure makes use 
of the results from the numerical calculations presented in Tenorio-Tagle et al. (1999), and 
illustrates in a simplified way the different elements that determine the properties of the Lya 
emission profile. We have identified 6 basic steps, as follows: 

i) Initially, when a star-forming episode starts, a central HII region begins to develop. At 
this phase, if the neutral gas surrounding the starburst region has HI column densities 
above 10^^~^^ cm~^, an absorption line centered at the systemic velocity of the galaxy 
will be visible, independently on the viewing angle. If the total HI column density 
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along the line of sight is higher than around 10^^ cm~^, a damped Lya absorption 
profile will be detectable. It is important to stress that during this early phase the 
Balmer lines will be strongest, due to the high ionizing flux produced by the most 
massive stars. 

ii) The situation changes drastically and becomes a strong function of viewing angle, once 
the mechanical energy released by the starburst is able to drive a shell of swept up 
matter to exceed the dimensions of the central disk. Then, upon the acceleration that 
this shell experiences as it enters the low density halo, it becomes Raylcigh-Taylor 
unstable and fragments. This event allows the hot gas (composed basically by matter 
recently processed by the starburst), to stream with its sound speed between fragments 
and follow the shock which now begins to form a new shell of swept-up halo matter. 
Another consequence of blowout is the fact that the ionizing photons from the recent 
starburst are now able to penetrate into the low density halo, and manage to produce 
an extended conical HIT region that reaches the outskirts of the galaxy. Given the low 
densities in the halo it is likely that this matter will remain ionized for a time that 
well exceeds the duration of the starburst activity. This also implies that some UV 
photons are at this stage able to stream freely into the inter galactic medium. The 
predicted expansion speed of this second shell formed in the halo would be around 
several hundred km s~^. 

An observer looking then at the starburst through the conical Hll region will be able 
to detect the strong Lya emission line produced by the central Hll region, centered at 
the systemic velocity of the galaxy. On the other hand, an observer looking outside 
the conical HIT region will still detect a broad absorption profile at any evolutionary 
state. 

iii) Sooner or later, recombination will begin to be frequent enough in the expanding shell. 
This will cause a strong depletion of the ionizing radiation which formerly was able 
to escape the galaxy after crossing the extended conical HII region. Recombination 
in the expanding shell will produce an additional Lya component, which the observer 
will detect blueshifted according to the expansion velocity of the shell, —Vexp- Clearly, 
if the shell were symmetrical and one could see the emission coming from the whole of 
it, a top-hat line profile would have to be added to the central HII region Lya emission 
line. If only a fraction of the shell is observed, we would expect to see 2 additional 
emission peaks centered at ±Vexp- 

iv) There are three efects that eventually lead to the trapping of the ionization front within 
the expanding shell. 
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— The first one is the increasingly larger amount of matter swept into the expanding 
shell, as this ploughs into the halo. 

— The growth of the shell dimensions also implies less UV photons impinging, per 
unit area and per unit time, at the inner edge of the shell. 

— Finally, in the case of a nearly instantaneous star burst, the production of UV 
photons starts to decrease drastically (as t~^) after the first 3.5 Myr of evolution. 

The trapping of the ionization front will yield the formation of a neutral layer at the 
external side of the expanding shell. All these effects lead then to an increasingly 
larger saturated absorption, as the external neutral layer will resonantly scatter the 
Lya photons. As discussed above, this absorption will be blueshifted with respect to 
the Lya photons emitted by the central HII region by —Vexp, leading to the formation 
of a P-Cygni profile with a fraction of the intrinsic Lya emission beeing absorbed. 

In addition, the profile will be contributed by the Lya radiation arising from the reced- 
ing section of the shell, both by recombination on the ionized layer, and by backscat- 
tering of the central Lya photons by the neutral layer. Since this contribution will be 
redshifted by 2xvexp with respect to the absorbing layer, it will be essentially free from 
resonant scattering and will be able to escape the region almost unaffected. 

Under such circunstances the resulting Lya emission line will be strongly depleted, 
showing a very asymmetric profile with a sharp blue edge and a line centroid shifted 
towards the red at a A different from the rest velocity of the host galaxy. The red wing 
of this profile will be furthermore broadened by the contribution from the receding 
shell. 

v) Under some conditions (for specific values of the shock velocity Vg and halo density 
Uhaio), the leading shock front may become radiative. An additional Lya component 
would thus arise as the shocked gas undergoes recombination. This should happen after 
a cooling time (U = l.5kT/ (inhaio-^); where T = 1.4 x W {Vs / lOOOkms-'^f , A is the 
coohng rate and k the Boltzmann constant) . Thus a combination of shock speed and 
density of the background halo may lead to a t\ smaller than the dynamical time and 
at that moment Lya emission will be produced as the shocked gas recombines. This 
additional Lya emission component would appear blueshifted at —Vg^p, and wouldn't 
be absorbed being ahead of the HI layers. 

vi) At a later phase, the shell will be completely recombined and have already substantially 
slowed down its expansion. Lya will show a broad absorption with only a small blue 
shift. The recombination lines from the HII region will be very weak. 
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The scenario here presented assumes a non-steady state of the star formation process. 
It imphes therefore either a (nearly) instantaneous starburst or the first phases of a star 
formation episode extended in time. In the case of a continuous star formation regime, a 
steady state where the birth and death of massive stars balances, so that the ionizing flux 
remains basically constant, is achieved only after the first 15 Myr of evolution (Cervino & 
Mas- Hesse 1994). If the process is active for much longer periods of time, the shell will finally 
blowout in the intergalactic medium, and no Lya emission lines with P-Cygni profile would 
be expected, but pure emission or absorption, depending on the orientation. 



5. General discussion 

Fig. 19 shows three GHRS low resolution examples corresponding to typical cases as 
described above. T12 14-277 shows a prominent and symmetric Lya emission line. The 
large equivalent width of the line (~ 90 A) points to a very young starburst. We postulate 
that this object would correspond to steps ii-iii above. In IRAS 0833-1-6517 the ionization 
front has already been trapped, so that a neutral layer has already developed, giving rise 
to a blueshifted absorption and the corresponding P-Cygni profile (step iv). Moreover, an 
additional blueshifted emission component can be easily identified in higher resolution images 
(see Fig. 9), as proposed in step v. Finally, SBS 0335-052 shows a very broad, damped Lya 
absorption line. This could hint to a very young object (step i), or to a geometrical effect, if 
we were observing from outside the ionization cone and through dense layers of static neutral 
gas. 

In this section we discuss in more detail the observed properties of the different objects 
considered, and how they fit within the proposed scenario. 

SBS 0335-052: SBS 0335-052 is a metal-deficient blue compact galaxy which is hosting a 
very strong star-forming episode. Its morphology and properties were described by Thuan, 
Izotov & Lipovctsky (1997). They identified 6 super star clusters containing an equivalent 
of ~4500 ionizing 07 stars. This galaxy is embedded in a large HI cloud with a mass of 
~ lO^M©. The HI column density along the GHRS aperture measured by these authors is 
very large, N(HI) = 7.0 x lO^^ cm'^. The observational properties of SBS 0335-052 would 
be consistent with step i: a very young starburst starting to ionize the surrounding medium 
while a significant amount of HI gas remains in front of the massive stars and is still static 
with respect to the central HII region. As a result, while a strong emission line spectrum 
is seen in the optical, the spectral region around Lya is completely absorbed, showing a 
damped absorption profile, similar to the ones computed in Fig. 14. 
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Kunth et al. (1998) identified the 01 and Sill absorption lines attributed to the neutral 

gas at the same redshift than the central HII region. Thuan & Izotov (1997) postulated 
the presence of two additional 01 and Sill absorption systems blueshifted by 500 and 1500 
km s~^, which might indicate that some layers of HI could be starting to be accelerated by 
the release of mechanical energy from the star burst, as expected in step i (see Fig. 18a). 

Fig. 20 shows the presence of a marginal excess emission in the wings of the Voigt profile, 
suggesting an intrinsically very strong Lya emission line. We have fitted the intrinsic Lya 
profile that would produce such an excess, after convolution with the absorption expected 
by the large amount of neutral gas, deriving a very strong emission with W{Lya) = 120 A, 
consistent with a very young and powerful starburst. 

IZwl8: The GHRS high resolution spectrum of IZwlS was first discussed by Kunth et al. 
(1994) and reanalyzed by Kunth et al. (1998). It ressembles very much the spectrum of 
SBS 0335-052, with a broad, damped Lya absorption profile corresponding to an HI column 
density N{HI) ~ 3.0 x 10^^ cm~^. IZwl8 is also a metal deficient blue compact galaxy 
surrounded by a dense cloud of neutral hydrogen. Kunth et al. (1998) showed that the 01 
and Sin interstellar absorption lines were at the same redshift than the central HE region, 
indicating that the neutral gas is unperturbed and mostly static with respect to the starburst 
region. 

As discussed above, no Lya emission is detected along the STIS slit, suggesting that a 
static neutral gas cloud completely covers the ionized region with a large column density along 
the fine of sight. This case seems in principle to be similar to SBS 0335-052. Nevertheless, 
while the starburst in this latter galaxy seems to be very young, there are hints for a continued 
massive star formation activity in lZwl8 over the last 10-15 Myr (Mas-Hesse & Kunth 1999), 
although at a rather weak rate. Moreover, large supershells and ionized filaments have been 
identified in this galaxy by Martin (1996). The observer might be looking in this case through 
the large amounts of static neutral gas expected outside the conical Hll region, which could 
be located along the supershells detected in the optical. 

T1214-277: The low resolution GHRS spectrum of T1214-277 has been described in detail by 
Thuan & Izotov (1997). This very low metallicity galaxy (~ Zq/ 23) hosts a very young and 
massive starburst. As these authors quote, this galaxy has a large Lya equivalent width (~70 
A). Campbell, Terlevich & Melnick (1986) measured also a very large H(3 equivalent width 
{W{HP) ~ 320 A within a 2" X 4" aperture). No any intrinsic extinction was measured by 
these authors from the Balmer lines ratios (but note that the Galactic extinction alone, with 
E(B — V)= 0.06, destroys around 45% of photons in the Lya region). The high equivalent 
widths and the absence of WR features (Pagcl et al. 1992) indicate that the starburst in 
T1214-277 has to be younger than 3 Myr, time at which the first massive stars enter the 
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WR phase at low metallicities (Mas-Hesse & Kunth 1999). 

These observational results indicate that T1214-277 could be at the end of step ii, as 
explained above: a conical HII region would have already developed, so that an observer 
looking through it could detect the strong Lya emission produced by the central HII region 
without significant distortion. 

In Fig. 21 we show a detail of the Lya emission profile in velocity scale. The profile 
is broadened, with 2 symmetric secondary peaks located at ~ ±300 km s~^ from the cen- 
tral emission peak. We postulate that the starburst in this galaxy may be indeed already 
entering step iii, so that recombination within the expanding shell produces the additional 
Lya components. Since only a fraction of the shell is seen through the GHRS aperture, the 
observer detects only 2 additional peaks centered at ±f e^^p- The velocity derived for the shell 
(300 km s~^) is within the expected range (few hundred km s~^) and is similar to velocities 
measured in the other starburst galaxies with P-Cygni profiles. 

Haro 2: Haro 2 is a well studied, rather metal-rich blue compact galaxy. As described by 
Lequeux et al. (1995) and Kunth et al. (1998) Haro 2 shows a prominent Lya emission 
line, with a clear P-Cygni profile. The metallic 01 and Sill interstellar absorption lines arc 
detected at around -200 km s~^ with respect to the central HII region, without any significant 
absorption feature detected at the redshift of the HII region. Haro 2 would be the prototype 
of our step iv, with a shell of neutral gas expanding at around 200 km s~^ with respect to 
the central HII region, ploughing into a low density, ionized halo without a screen of static 
neutral gas. 

Prom the STIS long slit spectral images (Pig. 1) we can identify all the features expected 
in step iv: 

• The emission decreases rapidly bluewards of the velocity at which the emission peak 
originated from the central HII region. The fiux goes to zero at basically the same 
velocity all along the slit. As discussed in Sect. 2, this is consistent with an expanding 
shell of neutral gas which has to be much more extended than the 8" over which the 
absorption is detected in order to behave as a plane-parallel slab. Prom the analysis 
of Ha long-slit spectroscopy Legrand et al. (1997) found evidence for the presence of 
an ionized expanding shell over a region around 25" in diameter (nearly 2.4 kpc). The 
lateral edges of the shell appeared at the same systemic velocity than the central HII 
region, decoupled, as expected, from the disk rotation. The Ha line is also clearly 
broadened at the edges, in good agreement with what would be expected from an 
expanding shell. 

As explained above, step iv foresees that as the starburst evolves, the trapping of the 
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ionization front will yield the formation of a neutral layer at the external side of the 
expanding shell, with no static neutral gas inside the shell. This would produce the 
kind of profile we detect on the central region of Haro 2 (see Figs. 3 and 15). 

• At the same time, recombination on the internal layers of the shell produces additional 
emission components at ±Vexp: in this case at ±200 km s~^. While the blueshifted 
Lya emission component produced in the internal layers will be destroyed by resonant 
scattering on the external neutral layers (both expanding at the same velocity), this 
blueshifted emission should be detectable in Ha. Legrand et al. (1997) concluded 
indeed that the Ha emission profile (when looking straight into the starburst region) 
showed evidence for two additional weak components at ±200 km s~^. 

These results indicate that a section of the shell receding at the same velocity should 
also be present. The internal part of the receding shell would produce Lya photons, as 
it generates Ha ones. And the neutral layers of the receding shell would backscatter 
in addition some Lya photons produced at the central HII region. Since the receding 
neutral gas lies at -Vexp with respect to the central HII gas it scatters the blue wing 
of the intrinsic Lya emission profile. A key point is that this receding shell lies at 
+2xvexp with respect to the approaching shell, so that the Lya photons emitted or 
backscattered there traverse unaffected the approaching neutral layers and reach the 
observer redshiftcd by +Vexp- As a result, we expect an extended, broad and weak 
component of Lya photons centered at around ±200 km s~^ in the case of Haro 2 (see 
the example in Fig. 17). 

This is exactly what we detect in our spectral images (Fig. 1). A weak extended 
emission is detected at both sides of the central Hll region on the slit, but especially to 
the upper side. It results from the Lya photons originating from or backscattered by 
the receding shell. We show in Fig. 4 the profiles of this emission at several arcseconds 
from the central HII region. The profiles show: 

— In both parts of the slit, SW and NE of the central HII region, the blue edge 
of the emission profile lies at the same velocity: the approaching neutral shell is 
destroying all photons blueshifted with respect to the velocity of the central HII 
region. 

— Broad cxtentions up to around ±500 km s~^. Profiles are more symmetric than 
the central Lya emission line and their peaks are clearly redwards of the Lya 
peak. Both profiles show quite extended red wings at basically the same veloc- 
ity. These wings originate by recombination and backscattering in the ionized 
and neutral layers of the receding shell, respectively. Since the scattering process 
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affects pfiotons witfiin a fiigfi range of velocities, as shown in Fig. 4, the backscat- 
tered component should cover basically the same velocity range measured on the 
absorption profile. We consider indeed that the good agreement between the ter- 
minal velocity of the P-Cygni profile (around —500 km s~^) and the terminal 
velocity of the extended Lya emission profile (around -|-500km s~^) supports our 
interpretation. 

IRAS 0833+6517: As discussed in Sect. 2, Lya emission in IRAS 0833+6517 has a more 
complex structure than in Haro 2 (see Fig. 7). In addition a secondary Lya emission is 
conspicuous on top of the damped absorption. Nevertheless, the overall properties of the 
Lya emission are very similar between both galaxies: 

• The flux vanishes bluewards of the same wavelength along the slit (hence the same 
velocity) whenever Lya emission is detected (10", 4.1 kpc). The extension of the 
neutral gas shell has to be very large in this galaxy, reaching a diameter close to 10 
kpc. 

• S and N of the central region where Lya is in emission, we detect again an extended 
broad and weak component redshifted with respect to the central HII region by around 
300 km s~^. The terminal velocity of the Lya proflle all along the slit is at around 
+700 km s"-*^. As in Haro 2, we associate this component to the emission originated at 
the receding shell, both by the ionized inner layers and by backscatter in the neutral 
external layers. 

We therefore argue that the superbubble phase in IRAS 0833+6517 corresponds to step 
iv in our model. The presence of an additional emission component on top of the absorption 
profile blueshifted by around —300 km s~^ with respect to the central HII region, supports 
our interpretation that an expanding shell is moving towards the observer at this velocity. As 
explained in step v, the leading shock in this galaxy should have become already radiative, 
undergoing recombination and originating the secondary Lya emission line we detect at the 
velocity of the shell. The fact that this secondary emission does not appear affected by 
neutral gas scattering indicates that the HI column density in front of the expanding shell 
(either static or also expanding) has to be rather small, as expected in our model. The 
smaller spatial extension of this secondary emission indicates that only a fraction of the 
whole surface of the leading shock is undergoing recombination. Conditions given in step v 
strongly depend on the density of the halo, which is most likely not uniform, explaining the 
patchiness of this emission. 
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We have overplotted in Fig. 22 our high resolution STIS spectrum of IRAS 0833+6517 
over the lower resolution GHRS one presented by Gonzalez Delgado et al. (1998). This 
comparison exemplifies how the analysis of P-Cygni profiles at low resolution may yield 
misleading results. The asymmetry of the intrinsic profile translates into an artificial redshift 
of the emission line peak when observed at low resolution (in addition to the effect discussed 
in Sect. 3, Fig. 16). This effect might be especially severe when analyzing galaxies at high 
redshift with medium or low resolution spectra. Furthermore, the lower resolution does not 
allow to identify the damped, black absorption plateau (which in this case is partially filled 
by the secondary emission peak), preventing an accurate derivation of the intervening HI 
column density. 

We want to remark that our model provides a simplified scenario for the formation of 
Lya P-Cygni profiles in starburst regions. The actual geometry and kinematical structure 
of some objects, hke IRAS 0833+6517, might be significantly more comphcated, although 
the overall view should still be valid. Let us highlight some observational properties that 
may be relevant: 

• Kunth et al. (1998) were not able to detect the interstellar 01 A1302.2 and Sill A1304.4 
lines in their high resolution GHRS spectrum. The analysis of the lower resolution 
data from Gonzalez Delgado et al. (1998) shows that both lines are blended within 
a very broad absorption profile, which in fact was extended over most of the Kunth 
et al. (1998) GHRS wavelength range. Gonzalez Delgado et al. (1998) report also 
very broad profiles on other interstellar lines, with FWHM up to around 1000 km s~^. 
They interpret these broad absorption profiles as evidence of large-scale motions of the 
interstellar gas around the starburst region. 

• The Silll A1206 absorption fine detected bluewards of 1230 A is split into two compo- 
nents (see Figs. 7 and 9). While the weakest component is centered at the redshift of 
the central HII region, the strongest one appears blueshifted by around —470 km s~^, 
which is higher than the expansion velocity attributed to the shell. Gonzalez Del- 
gado et al. (1998) also found from lower resolution GHRS spectra that the interstellar 
SiIIA1260 and Cll A1335 were blueshifted by 450 - 520 km s^^ with respect to the cen- 
tral HII region. We want to stress that no Lya emission at around —500 km s~^ can 
be detected. Therefore the high velocity gas has to be contained within the expanding 
shell. 

We conclude that while the bulk of the neutral gas shell seems to be expanding at 
around 300 km s~^, some gas might be moving within the shell in a rather chaotic way at 
significantly higher velocities. 
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Other galaxies: The model we have proposed could also explain the observational properties 
of the Lya line in other starbursts studied up to now. IIZw70 and Mrk36, discussed by 
Kunth et al. (1998), and Tol 65, presented by Thuan & Izotov (1997), would be similar to 
IZwl8 and SBS 0335-052: they show a damped broad absorption profile, and their interstellar 
absorption lines appear at the redshift of the central Hll region. The HI gas is therefore static 
along the line of sight with respect to the HII region, either because their starburst episodes 
have not been able to accelerate the gas (being very young or weakly energetic compared 
to the amount of neutral gas surrounding them) , or because we are observing them through 
the densest part of the HI disk. In this latter case Lya emission could possibly leak from 
other regions (for instance if an ionizing cone is allowing the photons to escape) and they 
might be scattered to our hne of sight. 

We are presently performing a survey with HST-ACS looking for regions in this kind 
of galaxies where Lya photons could be leaking. Preliminar results appear in Kunth et al. 
(2003). 

Other galaxies studied by Kunth et al. (1998), hke ESO 400-G043 and ESO 350-IG038, 
showing P-Cygni Lya absorption and blueshifted interstellar lines, would also be in the 
state described by our step iv, which seems to be the most frequent one among the galaxies 
showing Lya emission. 

We want finally to stress that the expanding layer causing the typical P-Cygni Lya line 
profile in starburst galaxies is extended and can span several kpc across a galaxy. As discussed 
above, in some objects this can be traced across several hundreds of pc away from the main 
burst of recent star formation. In the case of IRAS 0833 this spans more than 2.2 kpc across 
the galaxy. Our observations thus prove that this is not a well localized phenomenon, but 
rather a large scale one affecting the whole central region of these galaxies. Nevertheless, to 
explain the properties of this emission (lack of structure and dumpiness, velocity dispersion 
and intensity) the low density medium or halo of the host galaxy, into which the shell 
propagates, has to be very extended so as to allow for the existence of a continuous shell 
with similar properties (column density, velocity, etc..) over dimensions > 1 kpc. If this were 
not the case and the expansion had reached the galaxy edge, upon blowout the shell would 
have become Rayleigh- Taylor unstable and would have rapidly broken while destroying the 
absorbing layer. Thus, although the Lya asymmetric line profiles are indicative of strong 
outflows, these do not certify that supergalactic winds, venting the metals from the newly 
formed starbursts into the inter-galactic medium, have already developed. In fact if the 
expanding shell could reach the edge of its host galaxy, becoming Rayleigh- Taylor unstable 
and fragmenting, the Lya line profile to be observed would be the almost unattenuated line 
produced by the central HII region, without any absorption by neutral gas but not the kind 
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of profiles generally observed. 

6. Implications for galcixies at high redshift 

It is expected that all L* galaxies have undergone a phase of rapid star formation at a 
certain stage of their evolution. During this phase they should have become very powerful 
producers of ionizing photons. Very strong emission lines are therefore expected from these 
objects, and with rather large equivalent widths. For redshifts above 2; ~ 3, the Lya emission 
line would be observable in the optical range, making it therefore a priori an optimal tracer 
of star formation. Nevertheless, the first searches did not detect Lya emitters. It has been 
only recently that Lya emission from galaxies at redshifts above 2; ~ 3 has been detected, 
but at average luminosities weaker by 2 orders of magnitude than expected (sec Rhoads et 
al. (2000) for a complete set of references concerning Lya emission at high redshift). 

Our analysis of Lya emission line properties in nearby starburst galaxies has shown that 
there are several issues affecting the shape and intensity of the line, so that only in a small 
fraction of the cases the full intrinsic emission would be detectable. In this section we will 
extrapolate our results to high redshift galaxies, aiming to estimate the expected properties 
of Lya emission in these objects. 

Star-forming episodes detected in galaxies at redshift 2; ~ 3 seem to be significantly 
stronger than those detected in compact galaxies in the nearby Universe (see Shapley et 
al. (2001), Shapley et al. (2003) and references therein). Nevertheless, the properties 
of the Lya emission profiles are very similar to those observed in local starbursts. Lyman 
Break Galaxies exhibit a broad distribution of Lya strengths and profile-types, ranging from 
damped absorptions to pure emission, including also P-Cygni-like absorption and emission 
combinations (Shapley et al. 2003). Asymmetric P-Cygni profiles have also been detected 
by Frye, Broadhurst & Bcnitez (2002), Pettini et al. (2000, 2001), among other authors (see 
also Heckman (2001) and Rhoads et al. (2000)). 

Pettini et al. (2001) observing Lyman Break galaxies found velocity offsets of the in- 
terstellar absorption lines with respect to the velocity of the HII region, assumed to be at 
the systemic velocity of the galaxy. They found blue velocity offsets between approximately 
—200 and —400 km s^^, in three-quarters of the objects, with a median value of —300 
km s~^. More recently, Shapley et al. (2003) have obtained a composite spectrum by com- 
bining the data 811 individual Lyman Break galaxies. This spectrum reveals that the strong 
low-ionization interstellar features appear blueshifted with respect to the systemic velocity 
by an average of —150 ± 60 km s~^. Moreover, the mean velocity of the interstellar SilV 
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doublet is close to —180 km s~^, indeed very similar to the value measured for low ionization 
lines. As discussed above, the presence of blueshiftcd interstellar lines hints directly to the 
presence of neutral shells expanding at few hundred km s~^, as concluded by Pettini et al. 
(2001) and Shapley et al. (2003). On the other hand the simultaneous detection of SilV 
at roughly the same velocities suggests that the internal layers of the expanding shells are 
probably ionized, as discussed above. This range of velocities is consistent with the pre- 
dictions for a wind-blown shell, as we have seen. On the other hand, these authors found 
the Lya emission lines to be always redshifted with respect to the HIT regions. The offsets 
measured for Lya in the galaxies showing blueshifted interstellar absorptions span a range 
between around -1-200 and around -1-600 km s~^. The composite spectrum of Shapley et al. 
(2003) gives a mean rcdshift of -1-360 km s~"^ for the Lya peak with respect to the systemic 
velocity. As explained above this is exactly what one expects if the central starbursts were 
driving a neutral (or partially ionized) expanding shell (see Fig. 16). 

Pettini et al. (2000) show a high resolution, high signal to noise ratio spectrum of MS 
1512-cB58, a galaxy &t z = 2.73. They found a weak and asymmetric Lya emission line, 
redshifted by around +230 km s^^ with respect to the Hll region, on top of a damped 
absorption with A^(H1)= 7.5 x 10^*^ cm~^, blueshifted by —390 km s^^ (in agreement with 
the interstellar lines). At this rather high HI column densities the absorption profile is very 
broad, allowing only a small fraction of the intrinsic Lya emission line photons to escape. 
A careful analysis of their Fig. 4 shows the presence of a broad red wing, extending up to 
~ 800 km s~^ (velocity measured with respect to the HIT region). We postulate that this red 
wing component might be the emission originated at the receding part of the shell, either by 
recombination or by backscattering of Lya. photons from the central HII region, as discussed 
above. Wc expect this component to be broader the higher the HI column density, since 
backscattcr would affect photons at increasingly high velocities. The presence of the broad 
red wing provides additional support to the existence of a rather symmetrical expanding 
shell being energized by the central starburst. 

Asymmetric Lya emission profiles have also been detected by Frye, Broadhurst & 
Bcnitcz (2002) on a sample of eight lensed galaxies at high redshift, 3.7 < z < 5.2. Interstel- 
lar absorption lines have also been detected in some of these galaxies, but unfortunately the 
redshift of the central HE region is not available. They found an offset between the redshift 
of the interstellar hues and that of the Lya emission centroid in the range 300-800 km s~^. If 
we assume that the interstellar lines are tracing a neutral shell expanding at around 200-300 
km s~^ with respect to the central HII region, the Lya hues detected would be redshifted 
by up to around 400-500 km s~^ at most, in agreement with the predictions and with other 
galaxies. Note that Frye, Broadhurst & Benitez (2002) always quote the offset between Lya 
and the interstellar lines as to be the local redshift of the Lya profile, something which is 
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misleading since the shift with respect to the HII region has to be necessarily smaller. Fi- 
nally, most of the Lya profiles plotted by Frye, Broadhurst & Bemtez (2002) in their Fig. 13 
show a weak broadening of their red wing, in agreement with the model discussed in this 
work. 

Another important result obtained by Shapley et al. (2003) from their analysis of 
Lyman Break Galaxies is that the Lya emission strength increases as the kinematic offset 
between the Lya emission peak and the low-ionization interstellar absorption lines decreases. 
In addition, the Lya equivalent width in emission increases as the interstellar lines become 
weaker. We interpret these effects as an evidence that the lower the neutral gas column 
density, the smaller the relative redshift of the Lya emission peak and, correspondingly, the 
stronger the emission. This is in very good agreement with our above discussion (see for 
example Fig 16). 

We therefore conclude that the properties of the Lya emission profiles are very similar 
in local starbursts and in high redshift galaxies. Nevertheless, there are some differences 
between both kinds of objects, which have to be taken into account for a correct interpretation 
of the data: 

• Since the star-forming episodes observed at high redshift are very powerful ones, we 
expect intrinsically very strong Lya emission lines, with emission extending over several 
hundreds km s~^ at the base of the line. Even accounting for HI column densities 
around 10^° cm~^, photons at the red wing of the emission profile will be able to escape. 
Therefore, most of these high redshift galaxies experiencing strong starbursts should 
show Lya emission, albeit reduced by large factors from their intrinsic luminosity. 
Only galaxies with HI column densities above ~ 10^^ cm~^ will show a damped Lya 
absorption. 

• For the same reasons we expect the observed Lya emission lines to appear redshifted 
with respect to the systemic velocity of the galaxies by a broad range of velocities, 
up to several hundred km s~^. This effect should be smaller in local starbursts, since 
the Lya emission lines would be intrinsically weaker (and therefore narrower at their 
base). If no other line is available to measure the redshift, one should bear in mind 
that the derived value is only an upper limit of the galaxy redshift, and that the real 
one might be indeed smaller by several hundred km s~^. 

• The spectra presented by Frye, Broadhurst & Benftez (2002) show a feature which 
only appears for galaxies at high redshift: the Lya forest, that significantly lowers the 
average continuum bluewards of Lya. As a result, the blue wing of the absorption 
produced by the local expanding shell can not be identified in low resolution data. 
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On the other hand, the presence of asymmetric Lya profiles together with a disconti- 
nuity in the continuum (the average continuum redwards of Lya being stronger than 
bluewards of the hne) should allow to identify Lya in surveys aiming to estimate 
redshifts of galaxies at high redshift. 

• Stellar Lya absorption lines might be important at low resolution, artificially decreas- 
ing the observed intensity of the Lya emission line. At high resolution the stellar 
profiles can be detected and taken into account when measuring the line flux, but 
at low resolution both profiles might be blended and the emission line substantially 
hidden (see Valls-Gabaud (1993)). 

It is clear from the above discussion that using the Lya. emission line as a tracer of 
star formation episodes might lead to severe errors. We remark that both the measured line 
intensity, and thus the star formation rate derived from it, as well as its velocity, have to 
be considered as respectively lower and upper limits only. The problem might be especially 
significant in the case of star formation rates derived from Lya luminosities, which might be 
underestimated by more than an order of magnitude. 

From the above discussion we see that P-Cygni Lya profiles are predicted only when 
the supperbubble has entered the phase in which the ionisation front is trapped by the sector 
of the shell which is evolving into the extended halo. It has been noted in Tenorio-Tagle 
et al. (1999) that P-Cygni profiles are seen in galaxies that are on the higher luminosity 
side of the distribution {M < —18) in their sample. Similarly, at high redshift the Lyman 
Break Galaxies (hereafter LBG) exhibit similar profiles, as discussed above. Hence despite 
the much larger strength of the starbursts they are hosting, when compared to local ones, 
the HI gas is still present between the HII region and the observer. In particular this implies 
that in all these cases the ionizing radiation has to be trapped, and that a low density halo is 
surrounding these objects - as discussed above, without this halo the shell would have been 
disrupted and no P-Cygni profile would be produced. 

It is interesting to enlight this result with the constraints on the Lyman continuum 
(hereafter LyC) radiation from galaxies. Most of the previous observations of galaxies below 
the Lyman break (Leitherer et al. 1995; Hurwitz, Jelinsky & Dixon 1997) as well as recent 
measurements obtained from FUSE observations (Deharveng et al. 2001; Heckman et al. 
2001) show that the fraction of ionizing stellar photons that escape the ISM of each galaxy 
is small, if any. Most estimates led to an escape fraction < 5%. However these figures 
are in constrast with a LyC escape fraction larger than 50% reported by Steidel, Pettini & 
Adclberger (2001) obtained from a composite spectrum of galaxies ai < z > = 3.4. Indeed 
starburst galaxies at large redshift have a higher UV luminosity than most local starbursts 
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studied so far. This would imply more photons for ionizing the gas and probably an easier 
escape. However this difference remains a puzzle since one must reconcile two facts that seem 
to contradict each other. One is the leaking of the UV LyC and the other is the presence 
of P-Cygni Lya profiles. The very presence of these P-Cygni Lya profiles testifies that a 
substantial amount of HI is present, and hence should act as a screening agent against the 
escape of the LyC. One way out would be to assume that the HI layers are very perturbed and 
chaotic, so that they do not completely cover the massive stellar clusters. However the Lya 
emission that would escape from regions devoided of neutral gas should be overwhelming, 
resulting in a symmetric Lya emission profile. We predict therefore that LyC photons would 
be detected mostly on galaxies showing strong and quite symmetric Lya emission lines, 
rather than in objects showing typical P-Cygni profiles. As far as wc know, no relation has 
yet been established between the fraction of escaping LyC and the Lya emission profile in 
galaxies hosting starbursts at any redshift. 

There is a final point we want to remark. As discussed at the end of Sect. 4, the 
detection of Lya P-Cygni favours short-lived (or nearly instantaneous) bursts rather than 
steady state, long lasting star formation processes. After long-lasting star formation episodes 
the expanding shell would have finally blown-out in the intergalactic medium, and no P-Cygni 
profiles would be expected on the Lya lines. We would rather expect to detect pure emission 
or absorption, depending on the orientation of the line of sight with respect to the disk of 
the galaxy. In the case of short-lived starbursts, the usual concept of star formation rate 
(amount of gas mass being transformed into stars per unit time) becomes meaningless: in 
most starbursts, massive star formation has probably ceased already few Myr ago. It might 
be therefore misleading to derive the star formation rate density at high redshift from the 
analysis of galaxies which are forming stars in short bursts. To solve this inconsistency we 
propose to parameterize the strength of these star forming episodes by the total amount 
of gas transformed into stars, rather than by the rate of star formation, and to reevaluate 
accordingly the history of star formation in the early Universe. This analysis is out of the 
scope of this work, and will be presented in a forthcoming paper. 



7. Summary and conclusions 

Our findings have important implications for the study of Lya profiles in star-forming 
galaxies in particular, and for super-bubble evolution and galactic winds in general. We 
find that the P-Cygny profiles are extended, smooth and span several kiloparsecs covering a 
region much larger than the starburst and comparable or larger than the host galaxy itself. 
This strongly suggests the existence of an expanding super-shell generated by the mass and 
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energy loss of the star burst interacting with an extended low density gaseous halo. 

We identify six phases in the predicted evolution of a super-shell that correspond to the 
different types of Lya profiles observed. 

i) The early phase. In a very young starburst the recombination lines will be at their 
maximum strength and Lya will show a broad absorption centered at the systemic 
velocity of the galaxy. SBS 0335 could be in this stage. 

ii) The emission phase. The super-shell breaks, ionizing radiation escapes into the low 
density halo fully ionizing a "bi-conical" region. Lya is in emission at the systemic 
velocity of the galaxy. 

iii) Late emission phase. Recombination starts to dominate in the expanding shell pro- 
ducing an additional Lya component shifted by the shell expanding velocity. T1214 
would be at the end of phase ii, starting already this phase. 

iv) Shell recombination phase. Recombination dominates this phase characterized by a 
strong P-Cygni type profile in Lya. The absorption component is blue shifted by fg^rp) 
the expansion velocity of the shell. Lya emission is strongly depleted and the centroid 
of its profile is shifted toward the red. Backscattering and emission from the receding 
part of the shell may give raise to an extended red wing in the observed Lya profile. 

v) If the leading shock becomes radiative an additional unabsorbed and blue shifted com- 
ponent of Lya may be produced. 

vi) Late phase. The shell is almost completely recombined and has substantially slowed 
down its expansion. Lya will show a broad absorption with only a small blue shift. 
The recombination lines from the HII region will be very weak. 

These 6 phases and the associated 4 different profiles are able to describe, both qual- 
itatively and quantitatively, the variety of observed profiles in both the nearby HST 
sample and high redshift star-forming galaxies. Given that the probability of occur- 
rence of a given profile depends mainly on the age of the starburst (and also on the 
orientation) the observed fraction of systems in the different phases will provide im- 
portant information to critically test the predictions of this scenario. 

Our results stress the importance of the density and kinematical structure of the neutral 
gas surrounding an HII region on the detectability of the Lya emission. The observed 
properties of the Lya emission line will be a convolution of several effects. While the amount 
of dust alone determines only the absolute intensity of the emission line, but not affecting 
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its equivalent width, the kinematical configuration of the neutral gas is the driving factor for 
its final visibility, and its profile shape - from broad absorption to pure emission. The fact 
that several independent effects play such a significant role on the properties of the observed 
Lya emission lines explains the lack of correlations established in the past between the Lya 
emission strength and other properties of different starburst galaxies. The use of Lya as a 
tracer of star formation rate, and even as a redshift indicator for galaxies at high-z, should 
be done with care, always being aware of its limitations. 
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Fig. 1. — UV spectral images of Haro 2. Top: we have marked the position of the geocoronal 
Lya hne residuals, after background subtraction. The vertical line marks the expected 
position of the Lya wavelength at the redshift derived from the HII region. The pixels have 
been rebinned for better display. We show the angular scale along the spatial axis, as well as 
the corresponding spatial scale. The N and E arrows indicate the orientation of the sht on 
the sky. Bottom: Detail of the Lya, region. The image cuts have been selected to show the 
structure of the Lya emission. The displayed pixels correspond to physical detector pixels. 
The wavelength scale increases to the right on the X axis. 
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Fig. 2. — Spatial profile of the Lya line in Haro 2 on top of the UV continuum profile. Note 
the Lya emission extended to the SW (to the right in the plot), with almost no counterpart 
on the UV continuum. 
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Fig. 3. — Extracted spectrum of the central region of Haro 2 (corresponding to (y/4). The 
position of the geocoronal Lya line has been marked. Several interstellar absorption lines are 
detected bluewards of the Galactic Lya line. Note the broad but weak apparent absorption 
centered around 1240 A, of unknown origin. 
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Fig. 4. — Lya emission line spectral profiles in different regions of Haro 2. Strongest line: 
central region. Intermediate profile: extended emission integrated over a region of I'.'S cen- 
tered 2'.'3 SW of the nucleus. Weakest profile: extended emission over 0'.'7 centered at O'.'S 
NE of the nucleus. Note that the profiles red wing vanishes always at around 500 km s~^, 
independently of the strength of the line. Similarly, the blue wing goes to zero at around 
—200 km s~^ in the 3 regions. 
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Fig. 5. — Optical spectral images: Haro 2 (top) and IRAS 0833+6517 (bottom panel). The 
position of the HP and OIII hnes has been marked, as well as the angular and spatial scales 
along the Y axis. The pixels displayed correspond to physical detector pixels. Some velocity 
structure is apparent in the OIII lines on the Haro 2 image. The velocity shift between both 
peaks corresponds to around 700 km s~^. Such velocity offsets have not been detected on 
the Lya emission line. In the bottom panel we have marked with 'A' and 'N' the regions 
for which we have extracted the spectra shown in Fig. 11. Dark and white pixels are the 
residuals of regions affected by cosmic ray hits. 
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Fig. 6. — Optical spectrum of Haro 2, integrated over an aperture of &!3. Note the Wolf- 
Rayet feature at around 4686 A. 
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Fig. 7. — UV spectral images of IRAS 0833+6517. Top: we have marked the position of 
the geocoronal Lya hne residuals, after background subtraction. The vertical line marks the 
expected position of the Lya wavelength at the redshift derived from the HII region. The 
pixels have been rebinned for better display. We show the angular scale along the spatial 
axis, as well as the corresponding spatial scale. The N and E arrows indicate the orientation 
of the slit on the sky. Bottom: Detail of the Lya region. The image cuts have been selected 
to show the structure of the Lya emission. The displayed pixels correspond to physical 
detector pixels. The wavelength scale increases to the right on the X axis. The extracted 
spectrum is shown in Fig. 9. 
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Fig. 8. — Spatial profile of the Lya line in IRAS 0833+6517 on top of the UV continuum 
profile. The regions with the strongest Lya emission are not correlated with the location of 
the strongest UV continuum sources. 
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Fig. 9. — Extracted spectrum of the central region of IRAS 0833+6517 (corresponding to 
2'.'0). The position of the geocoronal Lya line has been marked, in the center of the broad 
Galactic absorption profile. We have also marked the position of the secondary emission 
peak, whose spatial distribution can be appreciated in Fig. 7. Note the splitting of the 
Silll A1206 absorption hue, with peaks at 1227 and 1229 A. Additional broad interstellar 
absorption lines are detected bluewards of the Galactic Lya line. 
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Fig. 10. — Lya emission line spectral profiles in different IRAS 0833+6517 regions. Strongest 
line: central region. Intermediate profile: extended emission integrated over a region of 2'.'3 
centered 3'.'0 N of the nucleus. Weakest profile: extended emission over 4'.'1 centered at 3'.'0 
S of the nucleus. Note again that the profiles extend over basically the same velocity range 
independently of their strength. The secondary emission profile, centered at —300 km s~^ 
appears only in the central region. 
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Fig. 11. — Extracted optical spectra of IRAS 0833+6517 for the regions marked in Fig. 5: 
A' region in the top panel, and 'N' region in the bottom one. 
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Fig. 12. — UV spectral image of IZwl8. Symbols as in Fig. 1. Note the broad damped 
absorption profile, blended with the Galactic absorption. Note also the lack of Lya photons 
in emission along the slit. 
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Fig. 13. — Extracted spectrum of the central region of IZwl8 (corresponding to 2'.'0). The 
abscissa axis is given in velocity scale to show the width of the absorption profile. Additional 
interstellar absorption lines are detected bluewards of Lya. 
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Fig. 14. — Voigtian absorption profiles computed for neutral hydrogen densities log(nHI) = 
18, 19.5, 20 and 21 cm~^ (from weaker to stronger absorption). We have assumed b = 20 
km s~^. 
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Fig. 15. — Expected Lya profiles. We have assumed an intrinsic Lya emission line origi- 
nating in the central HII region, at systemic velocity (V = km s~^). The plot shows the 
intrinsic emission profile (thick line, log(nHI) = km s~^) and the resulting profiles assuming 
a slab of neutral hydrogen moving at Vexp = —200 km s~^, with column densities log(nHI) 
= 19.5, 20 and 21 cm'^. 
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Fig. 16. — Detail of a resulting Lya profile showing that under certain circumstances the 
centroid of the observed Lya emission line can appear redshifted by several hundred km s~^. 
This redshift, nevertheless, is artificially originated by the absorption of the blue part of 
the profile, and should not be considered as a tracer of ionized gas outflows. Moreover, it 
is important to note that the width of the resulting line is much smaller than that of the 
intrinsic emission line. This example has been computed for log(nHI) = 20.3 cm~^ and Vexp 
= -300 km s-^ 
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Fig. 17. — Effect on the resulting profile of the additional Lya emission components. The 
first contribution originates at the inner part of the receding shell by photoionization and/or 
backscattering of Lya photons produced in the central HII region. The example has been 
computed for log(nHI) = 19.5 cm~^ and Vexp = 300 km s~^. The intensity of this contribu- 
tion has been assumed as 10% the intensity of the main Lya component originated in the 
HII region, as found by Legrand et al. (1997) for Haro 2. This component should appear 
redshifted at the velocity of the shell, in this case Vexp= +300 km s~^, and enhances the base 
of the red wing. The second contribution to Lya originates at the ionized region in front of 
the neutral expanding shell, appearing thus at the velocity of the approaching shell, Vexp = 
—300 km s~^. This component is not affected by neutral gas scattering, and should appear 
on top of the damped absorption profile. In the top panel we show the three contributions 
to Lya. The thick line is the emission line produced in the central HII region. The lower 
panel shows with thick line the resulting profile. The thin lines show the assumed absorp- 
tion profile and the convolved profile of the intrinsic Lya line alone. The effect of these two 
additional components on the total resulting profile is evident. 
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Fig. 18. — The basic model (see Tenorio-Tagle et al. (1999) for further details): evolution 
of the expanding shell generated by an HIT region and implications for the visibility of the 
Lya emission line, a) A massive starburst generates a central HIT region. The surrounding 
halo of neutral gas absorbs all photons with energy close to Lya ones, producing a damped 
absorption profile, b) At the beginning of the star-formation episode, the number of ionizing 
photons emitted by the central cluster of massive stars is very large, so that a fraction of them 
can escape the HII region and produce the ionization of the surrounding halo of the host 
galaxy. An observer looking straight through the ionization cone will detect a very strong 
Lya line centered at the rest velocity of the host galaxy. On the other hand, an observer 
looking at a higher angle will still see a damped absorption profile produced by the neutral gas 
on the galaxy disk, c) and d) The action of stellar winds and supernovae explosions generates 
an expanding shell that eventually will be able to undergo recombinations and emits also 
Lya photons. Furthermore, this large-scale expanding shell, driven by the mechanical energy 
released by the massive central starburst, is eventually able to trap the ionization front 
produced by the UV photons that escape the central HII region. This leads to the formation 
of two zones within the shell, expanding at the same speed. The inner zone is fully ionized 
and emits in Lya while the outer zone is neutral and capable of scattering and absorbing 
this radiation. The neutral gas in the approaching side of the shell leads to the formation 
of a Lya P-Cygni line profile, as discussed in the text. In addition, Lya photons produced 
by the central HII region and backscattered by the neutral layers of the receding shell will 
contribute with a low intensity, broad component redshifted by v^xp- 
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Fig. 19. — Observed Lya profiles illustrating the different cases discussed in the text: pure 
emission, with only weak absorption possibly produced by extinction; P-Cyg profile produced 
by a partially neutral expanding shell; and completely damped absorption produced by 
essentially static neutral gas. 
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Fig. 20. — Detail of the damped Lya absorption profile in SBS 0335-052. The thin line 
corresponds to the fitted Voigt absorption profile, including also the effect of the Galactic 
absorption. The thick line shows the convolution of the intrinsic emission line and the total 
absorption profile. It can be seen that the observed profile shows an excess with respect to 
the theoretical Voigt profile, detectable both on the red and the blue wings. We interpret 
this excess as the wings of the intrinsic Lya emission line generated in the HIT region, which 
become partially visible since the absorption at these wavelengths is not total, as evident 
from the Voigt profile. The HIT region Lya emission required to reproduce the observations 
would have an equivalent width around 120 A, as expected for a very young starburst. 
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Fig. 21. — Detail of the T12 14-277 Lya emission profile in velocity scale. The 2 peaks 
identified at ~ ±300 km s~^ from the central line peak could be attributed to the secondary 
emission from the approaching and receding parts of the young shell, according to steps ii - 
iii in the text. 
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Fig. 22. — Detail of the IRAS 0833+6517 Lya profile observed at two different resolutions. 
Note that low resolution spectroscopy can hide the presence of a blueshifted damped Lya 
absorption profile. 
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Table 1: Adopted properties of the observed HII Galaxies taken from the NASA Extragalactic 
Database, except for the metaUicity for which the references are indicated. 



Object 


RA(2000) 


Dec(2000) 


M(B) 


v(hel) 


Distance 


Scale 


12+log(0/H) 










km s~^ 


Mpc 


pc/" 




IRAS 0833+6517 


08 38 23.2 


65 07 15 


-20.8 


5730 


85.4 


413 


7.5^ 


IZwl8 


09 34 02.4 


55 14 32 


-14.0 


780 


10 


48.4 


7.22 


Haro 2 


10 32 31.9 


54 24 03 


-18.2 


1461 


19.5 


94.4 


8.43 



References. — 1- Margon et al. (1988); 2- Skillman & Kennicutt (1993); 3- Davidge (1989) 



Table 2. Journal of observations for the HST proposal 8302. The 52x0.5 slit was nscd in 
all cases. The UV spectra were obtained with the FUV-MAMA detector, at a central 
wavelength of 1222 A. The optical spectra were obtained with the CCD, at a central 
wavelength of 4300 A. The position angle of the Y axis is measured in degrees E of N. The 
coordinates indicated correspond to the center of the slit, as extracted from the STIS files 

headers (J2000). 



Object Obs. date Grating Integration time (s) Pes. angle (deg) Coordinates (RA, Dec - deg) 



IRAS 0833+6517 


2001-01-15 


G430L 


360 


171 


129.59666, 65.12083 


IRAS 0833+6517 


2001-01-15 


G140M 


1320, 3000, 3000 


171 


129.59666, 65.12083 


IZwl8 


2000-10-04 


G430L 


300 


-96 


143.50833, 55.24094 , 


IZwl8 


2000-10-04 


G140M 


1764, 3135, 3089, 3089 


-96 


143.50833, 55.24094 g 


Haro 2 7^1^(major axis) 


2000-02-21 


G430L 


300 


150 


158.13250, 54.40972 ' 


Haro 2 #1^ 


2000-02-21 


G140M 


1724, 3111, 3065 


150 


158.13250, 54.40972 


Haro 2 #2 (minor axis) 


2000-12-01 


G430L 


300 


-121 


158.13250, 54.40097 


Haro 2 #2 


2000-12-01 


G140M 


1724, 3111, 3065 


-121 


158.13250, 54.40097 



^Due to an error, the center of the slit was misplaced 31" to the N during this visit. 
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Table 3: Summary of observational results. The emission lines intensities have been measured 
through the same apertures in the UV and optical, so that their ratios reflect the intrinsic 
ratios in the gas. Lya has been measured extrapolating the UV continuum from longer 
wavelengths. The fluxes are given in units of erg s~^ cm~^. 



IRAS 0833+6517 Hare 2 



F[0//]3727 


6.0e-14 


1.2e-13 


F[C>777]4959 


2.3e-14 


6.9e-14 


F[O/77]5009 


6.5e-14 


2.06-13 


F{HP) 


2.0e-14 


7.7e-14 


W{H(3) (A) 


8 


75 


F{Lya) 


4.3e-14 


1.6e-13 


W{Lya) (A) 


6 


12 


F{sec.Lya) 


4.0e-15 




Aperture 


3'.'5x0'/5 


l'/5x0'.'5 



